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Abstract 

Mixed type (bright-dark) soliton solutions of the integrable N-coupled nonlinear Schrodinger 
(CNLS) equations with mixed signs of focusing and defocusing type nonlinearity coefficients are 
obtained by using Hirota's bilinearization method. Generally, for the mixed N-CNLS equations 
the bright and dark solitons can be split up in (N — 1) ways. By analysing the collision dynamics 
of these coupled bright and dark solitons systematically we point out that for N > 2, if the bright 
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amplitude dependent phase-shift and change in relative separation distance take place in the bright 
solitons during collision. However their counterparts, the dark solitons, undergo elastic collision 
but experience the same amplitude dependent phase-shift as that of bright solitons. Thus in the 
mixed CNLS system there co-exist shape changing collision of bright solitons and elastic collision 
of dark solitons with amplitude dependent phase-shift, thereby influencing each other mutually in 
an intricate way. 
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I. INTRODUCTION 



Solitons in coupled nonlinear Schrodinger (CNLS) equations have been the subject of in- 
tense studies due to their intriguing collision properties and their robustness against external 
perturbations. The study of physical and mathematical aspects of CNLS equations is of con- 
siderable current interest as these equations arise in diverse areas of science like nonlinear 



optics, optica 
physics 



jptical communic. 



communication, biophysics, Bose-Einstein condensates (BECs) and plasma 
In the context of nonlinear optics solitons arising as a result of 
competing focusing type nonlinearity and anomalous dispersion (diffraction) of pulse (beam) 
are called bright solitons [2I as they are well localized structures of light. The dark solitons 
resulting due to compensation between defocusing type nonlinearity and normal dispersion 
of pulse appear as localized intensity dips on a finite carrier wave background, and they are 
more robust than bright solitons The appearance of multicomponent CNLS type equa- 
tions as dynamical equations in various areas of physics and the subsequent studies 
on these systems lead to the identification of bright, dark, bright-dark, and dark-bright type 
solitons Even though there are a number of works on 



bright and dark soliton propagation and collision separately, still results are scarce for the 
study on bright-dark type soliton propagation and their collision dynamics. 

In the present paper we perform a study on bright-dark soliton solutions of mixed N-CNLS 
equations and examine the shape changing collisions of multicomponent bright solitons in 
the presence of dark components. We consider the following set of integrable mixed N-CNLS 
equations (in dimensionless form) 

N 

%> + gj,tt + 2(5^|g i | 2 )g J = 0, j = 1,2,..., N, (la) 
1=1 

where qj, j = 1, 2, . . . , N, is the complex amplitude of the jth component, the subscripts z 
and t denote the partial derivatives with respect to normalized distance and retarded time 
respectively, and the coefficients a^s define the sign of the nonlinearity. For convenience and 
without loss of generality, we define 07 's for this mixed case as 

o\ — 1 for I — 1,2, . . . ,m, 

= -1 for l = m + l,m + 2,...,N. (lb) 
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Recent theoretical and experimental studies show that the bright solitons in the focusing 
case (07 = 1, 1 = 1,2,..., N) undergo fascinating shape changing collisions characterized 
)y intensity redistribution, amplitude dependent phase-shift and relative separation distance 
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13], whereas the standard elastic collision of dark solitons occurs in the defocusing 
case (07 = —1, I = 1,2,...,N ) Later, it has been observed in Ref. [16] that 



bright-bright solitons of mixed CNLS equations ([I]) also undergo shape changing collision 
but of different nature. Very recently, it has also been shown that the bright solitons in the 
■nixed CNLS system exhibit periodic energv switching daring the shape changing collision 
process in the presence of linear couplings [171 ] . Now it is of further interest to examine how 
the bright solitons are influenced by dark solitons and vice- versa in this mixed CNLS system. 
The main aim of this paper is to investigate the nature of bright-dark soliton solutions of 
the mixed CNLS equations (1) and their collisions. 



The pioneering works of Makhankov et al. |18l in the context of Bose- Hubbard model and 

nn 

a few recent works jl9|, [20] on left-handed materials (LHMs) in nonlinear optics, suggest that 
the mixed CNLS system studied in our paper could be of considerable physical significance. 
Mixed 2-CNLS system can be obtained as the modified Hubbard model (Lindner-Fedyanin 
system) in the long-wavelength approximation by taking into account the electron-phonon 



interaction 



181 ] . Here the two component Bose-condensate is a mixture of two "gases" with 



attractive and repulsive interboson fields and the bright-dark soliton solution corresponds to 
the so called drop-bubble solution. A straightforward generalization of this mixed 2-CNLS 
system to an arbitrary number of fields given in Ref. [l8| . Exactly this is the system which 
we have considered in this paper (Eqs. (CEJ)) where the first m components have positive 
sign of nonlinearity and the remaining components possess negative sign of nonlinearity. 
The mixed 2-CNLS equations also arise as governing equations for an electromagnetic field 
propagation in LHMs with Kerr- type nonlinearity 19|. In fact, this case arises when the 
fields in the medium experience different types of nonlinearity, leading to o\ = —02 = 1 
or o"2 = — o"i = 1, corresponding to a medium of effective positive dielectric permittivity 
e e ff and effective negative magnetic permeability /i e // or vice- versa. This suggests that 
mixed N-CNLS equations could be a possible generalization of multiple electromagnetic fields 
propagating in LHMs with suitably chosen effective permittivity and effective permeability. 

Equations ([I]) with N = 2 can also be viewed as governing equation for two fields qi and 
q\ propagating in the anomalous and normal dispersion regimes, respectively (that is, the 
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self-phase modulation (SPM) coefficients are positive and cross-phase modulation (XPM) 
coefficients are negative in both the components). A possible physical realization of such 
type of nonlinearities is multi-field propagation in a quadratic medium with inefficient phase 



matching 



2l|. Another important physical realization of equations (1) arises in the context 



of Boson-Fermion gas mixtures. For example, the dynamics of two-species condensates is 
governed by mixed 2-CNLS equations for suitable choice of intraspecies (an, 022) and inter- 
species (012,0-21) scattering lengths [22|. These two-species condensates offer a wider range 
of possibilities, the main one being the possibility of having a negative interspecies scat- 
tering length. This possibility has been theoretically explored in the context of Feshbach 
resonance management and realized experimentally for boson-fermion mixtures 23j. Thus 
the interspecies interactions a 12 and a 2 ± can be tuned to be negative (attractive type) and 
positive (repulsive type) scattering lengths, respectively b y F eshbach resonance. It is also 
possible to choose the self interactions as a n > 0, a 22 < |24J. So studying N-CNLS equa- 
tions (1) of mixed type will provide a better understanding on the dynamics of multi-species 
condensates with suitably tuned scattering lengths. In addition, mixed 2-CNLS equations 
arise in BECs involving two isotopes of the same element, for example isotopes of rubidium 
(Rb 87 and Rb 85 ) 25|. In fact, multicomponent BECs support nonlinear waves which do 
not exist in single component BECs such as domain-wall solitons, bright-dark solitons, etc. 
Thus the mixed CNLS system which we have investigated could be of considerable physical 
relevance and significance in the context of nonlinear optics and matter waves. 

The plan of the paper is as follows. In Sec. II, we briefly present the Hirota's bilin- 
earization procedure for the mixed N-CNLS equations (1) to obtain exact soliton solutions. 
Section III is devoted to obtain exact one and two bright-dark soliton solutions for mixed 
2-CNLS, 3-CNLS and N-CNLS equations. The collision dynamics of bright-dark solitons in 
mixed CNLS equations is given in Sec. IV, where we have pointed out that if the bright 
solitons appear in more than one component then they undergo shape changing collisions 
characterized by intensity redistribution, amplitude dependent phase-shift and change in 
relative separation distance. We have also discussed the role of dark solitons on the shape 
changing collisions of bright solitons and also the effect of bright solitons on the dark soliton 
collisions. In Sec. V, we summarize the results of our study. Asymptotic analysis of mixed 
2-CNLS and mixed 3-CNLS equations are given in Appendices A and B, respectively. 
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II. BILINEARIZATION METHOD FOR INTEGRABLE MIXED N-CNLS EQUA- 
TIONS 



In this section, we briefly outline the procedure to obtain m-bright — n-dark soliton so- 
ution (m + n = N) of the mixed N-CNLS equations using Hirota's bilinearization method 



261 ] . We denote the soliton solution of Eqs. (1) in which the bright and dark solitons are 



split up in the first m components and the remaining (N — m) = n components, respec- 
tively, as "mixed soliton solutions" , for brevity. To start with, let us apply the bilinearizing 
transformation 

gU) 

Qj = -j-, 3 = 1,2, . . . ,m, (2a) 
h® 

%+m = — , I = 1,2, ...,n, (2b) 

to Eqs. (1), where g^''s and h^'s are arbitrary complex functions of z and t while / is a 
real function. Then the set of mixed N-CNLS equations given by Eqs. (1) reduces to the 
following set of bilinear equations: 

Di {9 {J) ■ f) = 0, j = l,2,...,m, (3a) 
D 1 (h®-f)=0, Z = l,2,...,n, (3b) 

tm n \ 

E^v^-E^^* . ( 3c ) 
3=1 1=1 J 

where Di = (iD z + D\ — A), D 2 = (D^ — X), * denotes complex conjugate and A is a constant 
to be determined. The Hirota's bilinear operators D z and D t are defined as 

^'■ b ^{r,-if)\m-wy^ t ^k'=^^ (4) 

Expanding gf^'s, h^'s and / formally as power series expansions in terms of a small arbitrary 
real parameter x, 

g U) = X g[ j) + X 3 g? + ..., j = l,2,...,m, (5a) 
h®= h$\l + X 2 h® + + ...), / = l,2,...,n, (5b) 
/= l + X 2 / 2 + X 4 /4 + ..., (5c) 



and solving the resultant set of equations recursively, we can obtain the explicit forms of 
g®% /i (/) 's and /. 



III. EXACT MIXED TYPE (BRIGHT-DARK) SOLITON SOLUTIONS OF MUL- 
TICOMPONENT MIXED CNLS EQUATIONS 



A significant feature of the system of integrable mixed CNLS equations (1) is that it 
admits a rich structure of soliton solutions like bright, dark, bright-dark soliton type solutions 
depending upon the boundary conditions. Here our focus is on bright-dark (or dark-bright) 
type solutions. In order to understand the nature of such mixed soliton solutions, their 
propagation and collision dynamics, we first obtain the solution for the N = 2 case. Then 
by extending the analysis to arbitrary N component case in the remaining sections, we point 
out that for N > 2, if the bright solitons are split up in two or more components, there 
occurs shape changing collision of bright solitons which is also influenced by the presence of 
dark soliton. 



A. Bright-dark soliton solutions of the mixed 2-CNLS equations 



a) Mixed (Bright-dark) one-soliton solution: 
In this case the bright soliton appears in the q\ component and the remaining component 
q2 comprises of dark soliton (or vice- versa). This case corresponds to the choice <j\ — 1 and 
a 2 = —I (or <Ti = — cr 2 = — 1) for which Eqs. (1) become 

m,z + q\,u + 2 (l<?i| 2 - l<?2| 2 )<?i = 0, 
iqi,z + q_i,tt + 2^\q 1 \ 2 - |g 2 | 2 )<?2 = 0. 
After restricting the power series expansion (jSJ) as 

- (7a) 



(6a) 
(6b) 



(i) 
X9i , 



(7b) 
(7c) 



and solving the resulting set of linear partial differential equations recursively, one can write 
down the mixed one-soliton solution explicitly as 
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Q2 



C\ e 



A x k m ^ ^sech 

k*+ib! 



R 



X _l_ e vi+vl+R 
6 



(8a) 
(8b) 



where 



771 = kit + i(k\ — X)z, ki = km + iku, = af R + ia^j , (8d) 

Ci = -{b\ + \)z + bit, A = 2|d| 2 , ci = ci R + icu, (8e) 

with the condition 

| Cl | 2 < \h -z&i| 2 . (8f) 

Here a^/c^and Ci are arbitrary complex parameters while bi is a real parameter. In the 
above equations and in the following sections, the suffixes R and / denote the real and 
imaginary parts, respectively. The above one-soliton solution is characterized by seven real 
parameters ct^, afj, km, ku, cm, cu and 61 along with the constraint (112). In the context 



of nonlinear optics the quantity Ai defined through the a-parameter can be viewed as 
the polarization vector of the light pulse/beam and Aikm as the amplitude of the bright 
soliton. By defining the quantities 6 and (j)^ as 6 = tan -1 (afj/ctm) and <pi = tan" 1 [(ku — 
bi)/km), respectively, Eqs. (8a) and (8b) can be rewritten as 

R 



qi = \l k\ R — IciPcos 2 ^^ e j6l sech 



k 1R (t-2k 1I z) . o 



x e ifci / t+j(fc2 Jj ,-fc2 / -2|ci| 2 )z^ ( 9a) 



q 2 = — Cie^ 1 ( cos^'j 1 '' tanh 



R 

km(t - 2k u z) + — 



+ zsin^ ) . (9b) 



Now the condition (l8fj) becomes k\ R > |ci| 2 cos 2 0^ 1 ' ) . It can be inferred from Eq. (l9al) that 
the intensity of the bright soliton increases with a decrease in the magnitude of the dark 
soliton parameter ci as shown in Fig. [TJ (all the quantities in this and rest of the figures 
are dimensionless) . This is a consequence of the particular type of cross phase modulation 
coupling given in Eq. ([6]). Note that the a-parameter influences only the central position of 
the bright and dark solitons and not the intensities. However, it can be inferred from Eq. 
(9a) that the a-parameters appear in the complex phase [e ld ] of the amplitude part, 
b) Mixed two-soliton solution: 

The mixed two-soliton solution can be obtained by terminating the power series expansion 
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as 



t(!) 



X9?+X 3 9 { s\ 



f = i + x 2 h + x'U 



(10a) 
(10b) 
(10c) 



After solving the resulting bilinear equations recursively, the explicit two-soliton solution is 
obtained as 



ft 

«2 



(11a) 



1 

D 



Q22 



+e 



(i) 



(lib) 



where 



and 



Y) — l _|_ g^l+^+'Rl _|_ g^l+J?5+5o _|_ e »?2+??J+<5(* _|_ e V2+ri2+R2 
_j_ e Vl +Vl +V2 +V2+R3 

rj j = k j t + i(k*-2\c 1 \ 2 )z, j = 1,2. 



(11c) 



(lid) 



In the above a[ , a<p , k\, k% and c\ are complex parameters and b\ is a real parameter (see 
below). Introducing the quantity 



(i) (1)* 
ol\ oq 



Hj (h + k*y 



Ci 



-1 



i,3 = 1,2, 



{k i -ib 1 ){k* j +ib 1 )_ 
the various other parameters in the expression (TTTT) are defined as follows: 

,<5o ,. „<5 f * 



e Rl = fi u , e R2 = n 22 , 



Jri _ (fc 2 - hY 



a 



(i> 



1 + 



= yUl2, e = /i 2 i, 
Icil 2 



o (1) 



(*2 - hf 



a 



(!)=» 



(/ci - ib 1 )(k 2 - 
2 



1 + 



c i 



(£4 - i6i)(fc 2 - ibx) 



(ki - «6i) 



and 



(A;* + i6i) 

R 3 _ I ^11^12/^21/^22 \ 



o (1) 



\ki - k 2 



1 + 



/^llA*21) 
/^22/^12; 

(fci - z6i)(fc 2 - ^1) 
(fcj +i6i)(^ +i6i) 

Icil 2 



^3 



(fei - i6i)(fc 2 - ibi) 



(12a) 

(12b) 
(12c) 

(12d) 
(12e) 

(12f) 



The above two-soliton solution (iTTj) is restricted by the conditions k 2 R +{kji—b\) 2 > \ci\ 2 , j = 
1,2, as in the case of one-soliton solution. The two-soliton solution (ITTj) is characterized by 
eleven real parameters a^, a^, ol£] , fcin, fci/, &2.R> &2/, Otis, and b\ . 



B. Bright-dark soliton solutions of the mixed 3-CNLS equations 

Let us consider the set of mixed 3-CNLS equations which corresponds to Eq. (TjQ) with 
N = 3. In its explicit form the set of mixed 3-CNLS equations reads as 

ki,z + + 2 ( a il<7i| 2 + c"2 k2 1 2 + 0"3k3| 2 )gi = 0, (13a) 
m,z + Q2,tt + 2^ai|gi| 2 + a 2 \q 2 \ 2 + o"3|<?3| 2 )<?2 = 0, (13b) 
m,z + Qz,u + 2^o"i|gi| 2 + cr 2 1 g 2 1 2 + cr^q^q^ = 0. (13c) 

In the above, the nonlinearity coefficients cr/s, j = 1,2,3 take the values either <j\ = a 2 = 
+1, and (T3 = —1 or G\ = +1 and cr 2 = cr 3 = —1. In Eqns. (13), in the context of BECs the 
components q±, q 2 and q% either denote the condensates of three isotopes of the same element 
(for example Rb isotopes) 2lj] or the hyperfine spin states of spinor BECs {27]]. This set of 
equations admits the following two distinct types of mixed soliton solutions. 

(i) 2-bright — 1-dark soliton solution: 

In this case the bright solitons are separated out into two of the three components and the 
dark soliton appears in the remaining component (with the choice, a% — a 2 = 1, and 03 = 
-!)• 

(ii) 1-bright — 2-dark soliton solution: 

In this type of solution the dark solitons appear in two components and the remaining 
component comprises of the bright soliton (with the choice, G\ — 1 and o 2 = a 3 = — 1). 
The procedure of obtaining these soliton solutions is similar to that of the mixed 2-CNLS 
equations. 



1. 2-bright — 1-dark soliton solution 

We present below the explicit forms of the obtained one- and two-soliton solutions, 
a) Mixed one-soliton solution: 

The mixed one-soliton solution, in which two bright solitons appear in the first two com- 
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ponents and the dark one appears in third component, is obtained by Hirota's method as 



I _|_ e v\+-nt+R ' 
A^i^'secli 



1,2, 



Ci e 



ViR + 2 



i ) e vi+vl+R 



I _|_ e vi+vt+R 



(14a) 
(14b) 

(14c) 



where 



,R 



A, 



E 



j=lV«l a 



(fci + fc!) 2 



I c il 



2k 



_R 

e 2 = 



|ci| 2 cos 2 ^»P 



1R 




a 



U) 



(1) 12 



Q: 



(2)m 



(14d) 
j = 1,2. (14e) 



Here the quantities rji and Ci are as defined in Eq. (8d) and (8e), respectively. A/s defined 
through the a- 7 and c\ parameters represent the polarization of the bright components. 
For the dark component, the parameter b\ denotes the direction of the background and 
Ci gives its amplitude. The one-soliton solution is characterized by nine real parameters 



a 



(1) (1) (2) (2) 



IRi a ll ' "lij> "1/ ' 



a; 



km, ku, Cm, cu , and b x and is restricted by the condition |ci| 2 < \k\ — 
ibi\ 2 . Now the role of a-parameters can be realized explicitly in the amplitude (intensity) 
of bright components and also through the non-trivial phase of all the components. This 
is shown in Fig. [2l In fact this has important consequences in the collision process as will 
be illustrated in the following sections. Thus the dark soliton part influences the bright 
part through the parameters c\ and b\ whereas the bright solitons influence the dark soliton 
phase (central position) through the a - parameters, 
b) Two-soliton solution: 

Following the Hirota's bilinearization method as in the case of iV = 2, here we obtain the 
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two-soliton solution as 



1 

D 



c^V 1 + a { 2 j) e m + e m+ ^ +V2+5 ^ + e ri2+ ^ +11 ' +5 ^ ) , j = 1,2, 



_|_ e ^l+''?r +»?2 +»75 +Q3 1) ^ 

e*« = (k 2 - k 1 )i2 11 /j2i(a { 2 j) X2i ~ Xn), 



,R 3 



\h - k 2 \ 2 nnfi 12 H2i^22{Xi2X2i ~ X11X22), J = 1, 2, 



where /ijj's are now redefined as 

1 



Xii 



(ki + k*)xu 
(hi + g) 



l c ll 



(fcj - ib\){k* + i&i 



z,Z= 1,2. 



(15a) 

(15b) 
(15c) 
(15d) 
(15e) 



(15f) 
(15g) 



The form of D is given as in Eq. (Illcj) and the expressions for e^", e^ 12 ' , e® 21 ' , e^ 22 
and e^s ' take the form as given in Eq. (12e) with the above redefinition of the /i«'s. 
a|jj, a£p, km, ka, i,j = I, 2, cm, cu and 61 are the fifteen real parameters which characterize 
the above solution. The nature of this two-soliton solution will be discussed in Sec. IVB. 



in 



2. 1-bright — 2-dark soliton solution 



Next we consider the case where the bright soliton appears in the qi component and the 
two dark solitons are found in the remaining two components (52,93)- This gives us the 
possibility of introducing two background fields cie^ 1 and c^e 1 ^ 2 ■ 
a) Mixed one-soliton solution: 



The corresponding one-soliton solution obtained by using Hirota's method is 



a 



I _|_ qVi+vI+R' 
A 1 k 1R e iVlI sech 



R 

ViR + g 



ci e 



Ql+i 



1 



ibi 



k* + ib t 



e m+vi+R 



I _|_ e m+vt+ R 



1,2, 



(16a) 
(16b) 

(16c) 
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where 



r ]l = k 1 t + i(k 2 1 - X)z, d = -(&? + X)z + b l t, A = 2(| Cl | 2 + |c 2 | 2 ), 1 = 1,2. 



e R = 



A 1 



(1)|2 



C 2 



\ki — ibA 2 \ki — ib 2 \ 2 



a 



(i) 



2k 



_R 

e 2 . 



1R 



(16d) 
(16e) 



Now the one-soliton solution is characterized by ten real parameters, a^,a^j, km, ha, 
cm, c u , c 2R , c 2I , bi and b 2 . 
b) Two-soliton solution: 

As in the previous section, here also we obtain the two-soliton solution as 

qi = 1 (aJV 1 + af'e* + e m+ ^ +r ' 2+Sl1 + e m+ ^ +111+5 ^ , (17a) 



1 



_|_g»72+T?|+Q22 _|_ g»71+>7l+»72+>72+Q 



Z = l,2, 



where 



(h - *i) 2 ^P, e fe = (k 2 - fc0 2 ^p. 



(i) 



(i) 

«2 



(17b) 



(17c) 



In Eqs. (17a) and (17b), the form of D, e Q ^ and e^' are same as in Eqs. (11c) and (12e) 
with 1 — 1,2. Also, the quantities e^ 3 , fin and p are redefined as 

e R3 = \ki - k 2 \ 



l«iv^ ' 



(fci + fc?) s 



1 rjt.. _ 7 



^ (fci - i&„)(A;,* + ib v ) 



, i,l = l,2 



(17d) 
(17e) 



and 



(17f) 



^ (fci - ^)(/c 2 - ib v ) 

The two-soliton solution is characterized by fourteen real parameters. Again we will study 
the nature of this solution in Sec. IVB. 



C. N-soliton solutions 



The above procedure of obtaining soliton solutions can be extended to three- and N- 
soliton solutions with some effort, though the analysis is cumbersome. In this work, we 
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restrict our analysis to the two soliton solution only as the N-soliton collisions represented 
by N-soliton solution in general take place pair-wise in soliton theory. Work is in progress 
in this direction and the results will be published separately. 



D. Bright-dark soliton solutions of mixed N-CNLS case 

After obtaining the two and three component mixed soliton solutions, the next natural 
step is to generalize the results to arbitrary N component case, where N = m + n. For this 
purpose, we consider the case where the bright solitons appear in the first m components 
and the dark solitons appear in the remaining n (= (N — m)) components. So the resulting 
mixed soliton solution can be denoted as m-bright— n-dark type soliton solution, as pointed 
out in Sec. II. 



1. m-bright— n-dark soliton solution 



a) One-soliton solution: 
The mixed one-soliton solution of the mixed N-CNLS case is found as 



where 



Here 



1 _|_ e vi+vt+R' 
Ajk 1R e irilI sech 



ViR + 



R 



c, e iCl 



Ql+r 



i _ / k-L-ibi \ Vl +ril+R 
1 kl+ih 1 



e R = 



A 3 = 



I _|_ e vi+vt+R 



j = 1,2,..., m, 
, 1 = 1,2, . 



1 - 



E \ci 
\h, - 



h - ik\ 2 



a 



2k 



e 2 . 



IR 



rjj = kjt + i{k) - \)z, d = -(bf + X)z + b t t, A = 2 ^ |q| 

i=i 

j = 1,2, ... ,m, and I — 1,2, ... ,n. 



(18a) 
(18b) 

(18c) 

(18d) 
(18e) 

(18f) 



The one-soliton solution is characterized by (2m + 3n + 2) number of real parameters, 
a iRi tt ?/i km, ku, cm, cu, bi, j = 1,2, ...,m, I = 1,2,..., n with the condition |q| 2 < 
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\ki — ibi\ 2 , 1 = 1, 2, . . . , n. 
b)Two-soliton solution: 

Generalization of the mixed two-soliton solution presented in the previous subsections for 
N = 2 and N = 3 cases yields the following m-bright— n-dark two-soliton solution of Eqs. 
(1) with arbitrary iV: 

1 



9iH 



1 

D 



ci 



j = 1,2, ...,m, 
e *C; _|_ e »7i+»7i+Qii _)_ e ii+»)2+Qi2 _|_ e V2+vt+Q 



(0 



Z = l,2, 



n. 



(19a) 



(19b) 



Here the denominator D is given by Eq. ( lllcl) . The quantities e Rl ,e R2 ,e s °, and e 5 o are as 
defined in Eq. (12b) but with the following redefinitions 



o (l) 



(M - ik) 



f^iji h j 1) 2, 



l(0 



(fci - z6 g )(fc 2 - gfci) 
(kl + ibi)(k^+ibi) 



(19c) 



e^ 3 = | - fc 2 | 2 /inA i i2A t 2i/i22[Xi2X2i - X11X22], 
1 



(19d) 



and 



Xip 



(k + k*)xi P 
(ki + k p ) 



i,p=l,2. 



(19e) 



The number of real parameters which characterize the two-soliton solution is (4m + 3n + 4). 



IV. SHAPE CHANGING COLLISIONS OF MIXED SOLITONS 



The fascinating property of the bright solitons of the integrable N-CNLS system with 
focusing type nonlinearity is that they exhibit shape changing collisions characterized by 
intensity redistribution, amplitude dependent phase-shift and relative separation distances 
{12I . I3I . 



16], which can then be used to construct collision based logic gates for optical 



computation 



28 



291 ] . In this and in the following sections we analyze such collision dynamics 



of bright solitons in the presence of dark solitons in the mixed 2-CNLS and 3-CNLS equations 
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and also its effect on the propagation and collision of dark solitons. In this regard, we perform 
an asymptotic analysis of the two-soliton solutions for mixed 2-CNLS and 3-CNLS equations. 

A. Asymptotic analysis of mixed two-soliton solution of mixed 2-CNLS equations 

To start with, we consider the collision properties associated with the mixed two-soliton 
solution (11) of the mixed 2-CNLS equations Following this we carry out the analysis 
for mixed 3-CNLS equations. Without loss of generality, we take kjR > and ku > k 2 i, 
kj = kjR + ikji,j = 1,2, and obtain the asymptotic forms of two colliding solitons (say Si 
and S 2 ). Similar analysis can be carried out for other choices of kjR and kji also. Using the 
expression (11) for the bright-dark two soliton solution of the mixed 2-CNLS system ([6]), we 
carry out a detailed asymptotic analysis in Appendix A for the two soliton collision process. 
Based on this analysis we identify the following: 

1. The role of dark soliton on bright soliton collision: 

The amplitudes of the two solitons S\ and S 2 before (after) interaction are given by 
A\~km (Al + ki R ) and Al~k 2R (Al + k 2R ), respectively, in the qi component. Forms of A 3 ^,j = 
1,2, are given in Appendix A (see Eqs. (Ale), (A2c), (A3c) and (A4c)). By rewriting these 
forms one can show that the intensities of solitons before and after interaction are same 
(elastic), (i.e.) \A{~~\ = j = 1,2, even though the complex amplitudes differ in phase. 

Also, the appearance of the background parameter l ci in the expression for the bright soliton 
amplitudes before and after collision (see Eqs. (Ale), (A2c), (A3c) and (A4c)]) shows that 
this parameter influences the bright soliton amplitudes throughout the collision. 
Typical bright-dark soliton collision in the mixed 2-CNLS system is shown in Fig. [3] for the 
parametric choices kx = 1 + i, k 2 = 2 — i, \c\\ = 0.56, bi = 0.2, = l,ot 2 = 1 + i. For 
better understanding, first we plot the intensity profiles showing the bright soliton elastic 
collision scenario in the absence of dark soliton parameter (c± — 0), at z — —4 and z — 4 . 
This is shown in Fig. H](a). Then the same bright soliton collision is plotted for the above 
parametric choices as in Fig. [3] at z = —4 and z = 4, that is in the presence of the dark 
component, in Fig. H](b). From these two figures, we observe that due to the presence of 
dark soliton the amplitudes of the colliding bright solitons are reduced by the same amount 
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throughout the collision. 

The two colliding solitons Si and S 2 suffer phase-shifts $1 and $2, respectively in both the 
bright and dark components. These phase-shifts for the bright and dark solitons are given 
by the expression 

i? 3 — R 2 — R\ 



$1 = — $2 




(20a) 



(k l + k*y(k 2 + kiy [1-^] 

where P\ = k\ — ib\ and P 2 = k 2 — ib\. Note that the phase-shifts appearing here (Eq. (20a)) 
and in the following (see Eqs. (23a) and (26)) are real quantities as the terms appearing 
in the argument of log function are products of complex conjugates. The role of dark 
component comes into picture through this phase-shift due to the explicit appearance of the 
background parameters (ci, &i). Notice that in the absence of the dark component {c\ = 0), 
the phase-shift $1 reduces to the standard phase-shift experienced by colliding solitons in 
scalar nonlinear Schrodinger equations (NLS). Also, it is important to notice that the a- 
parameters have no effect on the phase-shift in the mixed 2-CNLS case. However, this is not 
true in the case of mixed N-CNLS case for more than two components (N > 2) as will be 
shown in the following sections. The relative separation distances between the two colliding 
solitons before interaction 

_ RzkiR ~ Ri(k m + k 2R ) 

t 12 — r>; ; (ZUDJ 

2ki R k 2R 



and after interaction 



+ _ R2(k 1R + k 2R ) - R 3 k 2 R , v 

12 ~ 2k 1R k 2R [ C) 



also remain unaffected by a-parameters. 



2. The role of bright soliton on dark soliton collision: 

The asymptotic expressions [Eqs. (Alb), (A2b), (A3b) and (A4b)] reveal the fact that 
the collision among dark solitons is elastic as the intensities of colliding solitons remain the 
same before and after interaction. From Fig. [3] it can be observed that the bright component 
parameter 'a' has no influence at all either in the amplitude or in the phase-shift of the dark 
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soliton during collision. Thus the inclusion of a-parameters in the bright soliton solution 
(see Eq. (11)) does not affect dark soliton collisions. 

B. Asymptotic analysis of mixed two-soliton solution of mixed 3-CNLS equations 

The next natural step is to study the collision process in the mixed 3-CNLS equations and 
one can generalize the results to the N-CNLS equations with arbitrary N. The asymptotic 
expressions of the solitons corresponding to the N = 3 case are presented in Appendix 
B. This 3-CNLS system admits two distinct types of solution as mentioned in Sec. IIIB. 
First let us consider the 2-bright— 1-dark soliton collision and the 1-bright— 2-dark soliton 
collision. 



1. 2-Bright— 1-Dark Soliton Collision: 



We analyze two important physical quantities, namely (i) intensity and (ii) phase-shift 
of both bright and dark solitons. 
(i) Intensities of bright and dark solitons: 

Analysing Eqs. (Blc), (B2c), (B3c) and (B4c), we find that the amplitudes (intensities) 
of bright solitons before and after interaction are different. In fact, the intensities of the 
colliding bright solitons before and after interactions can be related through the expression 



\A?\ 2 =\Ttf\A l f\\ j J = 1,2, 



where the superscripts /± represent the solitons designated as Si and S 2 , a.t z 
transition intensities are identified from the Appendix B as 

\l-K 2 (a?/a?)\ 2 



IT 



IT 



1|2 
/ I 

2 12 
i I 

Xn 



11 — 

II - K\K 2 \ 



\\-K X {af */a?)\ 2 

Xn X22 
, ^2 = , 

X21 Xl2 

(hi + kt) 



1,2, 



c i 



, i,l = 1,2. 



(21) 
±oo. The 

(22a) 
(22b) 
(22c) 
(22d) 



Note that the transition amplitudes T-'s now are also functions of dark soliton parameters 



Ci and &x- On the other hand, from Eqs. fIBlbl) . flB2bl) . flB3bj) and flB4bj) we find that the 
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intensities of dark solitons remain unchanged due to the collision process, 
(ii) Phase-shift of bright and dark solitons: 

The amplitude dependent phase-shift $i( = i?3 ~ i ^ 2 ~* Rl ) for soliton Si can be expressed in 
the present case as 



$! = 1 log 



where 



l*i - h\ 2 (, u ( (fa + fa)(fa + fa*) [i - f^] [i - fflf] 



(fa + A3) (fa + A?) I Ufa + fa) (fa + fa) [i - ^] [1 - gjr] 



(23a) 



/ (1) (1)* . (2) (2)*w (1) (1)* . (2) (2)*s 
jj = K a 2 + «1 <*2 )(«2 <*1 +«2 «1 ) / 23b ^ 

(|af ) | 2 + |af ) | 2 )(|4 1) | 2 + |4 2) | 2 ) 

Note that the soliton S 2 experiences an exactly opposite phase-shift $2 (= — 3>i). The above 
phase-shift in turn results in the following change in relative separation distance between 
the solitons before and after collision, 

At 12 =t- l2 -t+ 2 = ^\ hR h l . (24) 

The expression for $1 is given in Eq. (I23al) clearly indicates that now the phase-shift depends 
on a-parameters as well as c\. Thus the phase-shift and ultimately the relative separation 
distance between the solitons can be altered during a two-soliton collision process for a given 
combination of fa, C\ and b\, by just varying the a parameters and as a whole the combined 
soliton profile gets altered. 

This kind of collision scenario is shown in Fig. [5] The corresponding intensity plots show 
that the bright solitons undergo shape changing collisions characterized by intensity redistri- 
bution, amplitude dependent phase-shift and relative separation distance for the parametric 
choice fa = 1 + i,k 2 = 2 - i,b 1 = 0.2, |ci| = 0.56, af } = l,af ] = (32 + z80)/89, a} ] = 
1, a 2 = 1. The solitons (say S± and S 2 ) are well separated before and after collision in both 
the components q% and q 2 . In the q\ component the intensity of soliton S\ gets suppressed 



while that of soliton S 2 is enhanced after interaction, whereas in the q 2 component it gets 

(1) (2) 

~w / (2) 



a (l) a (2) 

reversed. Such shape changing collisions occur for -fey 7^ -fey, which is quite general. But 



when we choose -fey = -fey, the two solitons exhibit elastic collision. It is instructive to note 



a 2 a 2 

that although the dark solitons appear in the q% component, they indirectly influence the 
shape changing collisions through the carrier wave background parameters C\ and hi. From 
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the asymptotic analysis it follows that these background parameters influence the intensities 
of the colliding bright solitons before and after collision by different amounts through their 
explicit appearance in the transition intensities (see Eq. (1221) ). However the nature of the 
collision is unaltered. We present below a detailed discussion to get a clear picture about 
the influence of dark solitons on bright soliton collision and vice-versa. 

(a) Effect of dark soliton on the intensity of bright soliton: 

For a better understanding, we present the shape changing collision of bright solitons (i) in 
the absence of dark component, that is C\ = 0, in Fig. [6] and (ii) in the presence of dark 
component, that is C\ ^ 0, in Fig. [7]for the above mentioned parametric choices. From these 
figures we observe that, in the presence of dark component, the intensities of solitons Si in 
the q\ component is decreased (increased) before(after) collision (as compared with Fig. [HD, 
but not by the same amount. Thus the effect of the dark component on the intensity of 
the bright soliton before its collision is different from the effect on the intensity after the 
collision. But in the q2 component there occurs decrement of intensity in both the solitons 
before and after collision due to the presence of dark soliton. This confirms the fact that the 
presence of dark solitons indeed influences the cross phase coupling between the two com- 
ponents qi and q2 which in turn affects the energy redistribution between those components 
as observed from figures. 

However, the nature of collision, that is enhancement (suppression) in Si and suppression 
(enhancement) in S2 in q\{q2) component during shape changing collision process, is still 
preserved. We also notice that in the case of standard elastic collision process resulting for 

a W (2) 

the choice -hr = ^W, the role of the c\ parameter on the amplitudes of colliding bright 
solitons before and after collision is same. This is shown in Fig. [8] and Fig. [9j for c\ = 
and c\ 7^ respectively. From the figures we observe that the intensities of bright solitons 
before and after collision are affected by the dark component by the same amount. 

(b) Collision behaviour of dark solitons in the presence of bright solitons: 

The intensities of dark solitons in the qs component are unaffected during the collision in the 
presence of the bright solitons in the q± and q2 components. This is obvious from the analytic 
expressions (IBlbl) . flB2bl) . flB3bj) and flB4bj) . This type of collision scenario is shown in the 
third figure of Fig. [5] (see also Fig. [7]). The analysis presented in Appendix B reveals the 
fact that the a parameters do affect the interaction of dark solitons through the shift in cen- 
tral position of the solitons which ultimately changes the separation distance between them 
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after collision. This phase-shift and the resulting change in the relative separation distance 
between the solitons can be obtained from Eq. (123aj) and Eq. (1241) . respectively. In fact, 
the change in relative separation distance becomes more significant and displays interesting 
propagation and collision dynamics of solitons when the soliton velocities are moderately 
different and kjRS are equal. For illustrative purpose, we consider the propagation of such 
composite two dark solitons arising for the choice k\ = 0.6 — i, k 2 = 0.6 — 0.5i, b\ = 0.2, |ci| = 
0.56, a{ 1} = 1, af ] = i, = (22/55) - 45i, a® = 1 in Fig. QUI at z = -5 and z = 5. The 

/ a (l) Q (2) \ 

analytic expression corresponding to the above choice ( k 1R = k 2R , -jtj ^ -^y ) is given by 



e%cosh(A + m + e^^cosh^ + SszOn) + e ^ 21 cos(C + Q^i^n 



Qll+Q'22 



Q12+Q21 



where 



e^cosh(A + e 5 4 52 cosh(5 + ^f^) + e^costC 



2 



(25a) 



.4 
B 
C 



2k 1R [t - z(k u + k 2I )}, 

2k 1R [z(k u - k 2 i)), 

(k 2I - k XI )t + (klj - k 2 2I )z. 



(25b) 
(25c) 
(25d) 



ad, 



3JJ and all the parameters in 



Such a soliton solution can be viewed as soliton complex 
Eq. (125aj) are defined in Eq. (15). Fig. [TUlfa) shows that the dark soliton complex varies 
its profile during propagation within a finite distance. This is a consequence of change 
in the relative separation distance due to the presence of a-parameters. However for the 

(1) (2) 

same choice of fe,'s, 61, and c\ but with = 5!i the dark soliton complex does not 
vary its profile after propagation through the same distance. This is shown in Fig. ITUlfb) . 
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31] . The bright 



This kind of behaviour is a striking feature of multisoliton complexes 
counterparts of this dark soliton also form bright soliton complexes of variable shape and 
such bri ght soliton complexes in CNLS equations with focusing nonlinearity are discussed 



in Refs. 



m 



30 



31 



32] . Work is in progress on such bright-dark multisoliton complexes. 
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2. 1- Bright— 2- Dark Soliton Collision: 



Two soliton solution given by Eq. (17) describes the collision of solitons in which bright 
soliton collision takes place in the qi component and the dark soliton collision occurs in 
the q 2 and g 3 components. The asymptotic expressions for the colliding solitons before and 
after collision are given in Appendix B. Equations (B5), (B6), (B7) and (B8) show that 
the intensities of bright solitons before and after collisions are same , (i.e.) \A l f\ = \A l ^~\, 
I = 1,2. Similarly the intensities of the dark solitons are also unaltered after collision. This 
indicates that the bright as well as the dark solitons undergo elastic collision. However there 
occurs a dark soliton parameter dependent phase-shift due to collision. The phase-shift of 
soliton Si in both the bright and dark components is given by 



( 


ki - k 2 


2 


1 1 Y^ 2 KI 2 

L Z^v=l (ki-ib v )(k 2 -ib v ) 


2 \ 


(kl + k*) 2 (k 2 + kl) 2 


1 \c v \ 2 

Z^v=l (k 1 ~ib v )(k 2 -ib v ) 


2 

/ 



(26) 



Similarly soliton S 2 undergoes a phase-shift $ 2 = — The dark component parameters 
Ci,c 2 ,bi and b 2 influence the phase-shift whereas the bright component parameter a does 
not alter the phase- shift. 

The above asymptotic analysis of two soliton solution of 2-CNLS and 3-CNLS can be 
extended straight forwardly to N-CNLS equations, with arbitrary N. By generalizing the 
above study we also point out that for mixed N-CNLS equations (1), the bright and dark 
solitons can be split up in (N — 1) ways, starting from (N-l)-bright— 1-dark soliton solution, 
(N-2)-bright— 2-dark soliton solution, up to 1-bright— (N-l)-dark soliton solution. We also 
arrive at an important conclusion that for m-bright— n-dark two-soliton solution, where 
m + n = N, the shape changing collision of bright solitons takes place only when m > 2. 



V. CONCLUSION 



To conclude, we have obtained the explicit mixed type (bright-dark) soliton solutions 
for the multicomponent mixed coupled nonlinear Schrodinger equations using Hirota's bilin- 
earization method. In particular, we have shown that these coupled bright and dark solitons 
possess rich structure and become more general than individual bright /dark solitons. Next, 
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our analysis of their collision behaviour reveals the fact that there exist elastic collisions of 
bright as well as dark solitons for the two component case. We observe that for this case the 
bright solitons are affected uniformly by the dark solitons. The bright soliton parameters 
a's have no effect on the phase-shift whereas it is influenced by the dark soliton parameters 
C\ and b\. The important observation of this study is that for more than two components, 
if the bright solitons appear in at least two components, then those bright solitons undergo 
shape changing collisions characterized by intensity redistribution, amplitude dependent 
phase-shift and relative separation distances, but their counterpart (dark solitons) exhibits 
only elastic collisions but with amplitude dependent phase-shift. This identification can 
find potential applications in optical as well as matter wave switching devices where the 
switching is performed through shape changing collision of solitons. Further, in contrast to 
the N component Manakov system here the phase-shift of bright as well as dark solitons is 
characterized by dark soliton parameters c\ and b\ in addition to the a parameters. Also, 
we observe that the dark soliton parameter c\ influences the intensity of bright solitons by 
different amount when the bright soliton parameters a's are such that -jw 4 -757 while their 

(!) (2) 

amplitudes are affected by dark solitons by the same amount, for = ^W. One more no- 
ticeable observation is that the dark solitons vary their profiles depending on a parameters 
during propagation as in multisoliton complexes 
the study will give further insight into the bright-dark paired solitons, soliton complexes 
formation, collision in Boson-Fermion mixtures and in nonlinear left-handed materials and 
their applications in switching devices. 
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32| . The various results obtained from 
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APPENDIX A: ASYMPTOTIC ANALYSIS OF BRIGHT-DARK TWO- 
SOLITON SOLUTION OF MIXED 2-CNLS EQUATIONS 



In the limit z — > ±00 the two-soliton solution (11) takes the following asymptotic forms 
for the choice k 1R , k 2R > and k u > k 2 i- For other choices of fcj's, i = 1,2, similar analysis 
can be made. 

a) Before collision (limit z — > —00): 
(i) Soliton 1 (r) 1R ~ 0,r] 2R -> -00): 



where 



g 2 ~ - Cl ^Ci+tf') 



cos^^tanh ^771^ + ) /sino. 



(1)* / i., „ 



k\ R — |C!| 2 COS 2 0^, 



<t>\ 



(i) 



tan 



if hi- h 



"1R 



(1) 



(Ala) 
(Alb) 

(Ale) 
(Aid) 



(ii) Soliton 2 (rj 2R ~ 0,i] m -> 00): 



9i - ^ A; 2 ; J e M?2 'sechU 2i? + 



(i* 3 - #i) 



g 2 — ci e 



(fcj + i&i) 



COSi 



tanh ^?7 2R + 



(i?3 — 



where 



and 



2- 



Q. 



tan" 



if hi- h 

h R 



In the above expressions, 



+ isin02 1 ' ) 



(A2a) 
, (A2b) 

(A2c) 

(A2d) 



X = (£4 - fc 2 )(A;i + k* 2 )(kl + i6i)([(fci - i6i)(A; 2 * + ih) - \ Cl \ 2 ][(h - ih)(h - ih) + | Cl | 2 ])2. 
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b) After collision (limit z — > oo): 
(i) Soliton 1 (rj 1R ~ 0,r] 2R -»■ oo): 

, (Cl+ ^,W) (fc 2 - L^(D, 



g 2 — Ci e 



(AS + 



where 



1+ 



1 



in, u f . (R3 — R2) 
cos(p\ 'tanhirjm H 



Y 



Here 



+ isin0^ 



(i) 



(A3a) 



, (A3b) 



(A3c) 



Y = (ki - k 2 )(k 2 + kt)(k* + ibi)([(ki - ibjfa - i&i) + IciHpJ + ibjfa - ih) - \ Cl \ 2 ]^ 



(ii) Soliton 2 (r? 2J? ~ 0,r] 1R -> -00): 



where 



?2 



A 2+ 



-ci e 



i(Ci+4 1) ) 



cos^^tanh ( 77^ + ^ ) + zsi 



tsm< 



1 



J~^\/ k 2R - |ci| 2 cos 2 ^ (!l 



(A4a) 
(A4b) 

(A4c) 



Note that in the above expressions though A l + ^A l ~, j, I = 1, 2, |Aj+| = |Aj _ | an d thereby 
confirming that no intensity redistribution occurs in the case of bright soliton and also for 
dark soliton (see Eqs. (Alb), (A2b), (A3b) and (A4b)). The phase-shift $1 ($ 2 ) of soliton 
Si (S 2 ) in the bright and dark component obtained from the above asymptotic expressions 
is given by Eq. (20a). 



APPENDIX B: ASYMPTOTIC ANALYSIS OF BRIGHT-DARK TWO-SOLITON 
SOLUTION OF MIXED 3-CNLS EQUATIONS 

1. 2-bright — 1-dark soliton solution 

Here also we assume k± R , k 2R > and ku > k 2 i, then the two soliton (2-bright— 1-dark) 
solution (15) takes the following forms asymptotically (z — > ±00). 
a) Before collision (limit z — > —00): 
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(i) Soliton 1 (r) 1R ~ 0,r) 2R -> -oo): 



/.•|/;s<x-l, ( //,/, • ^ ) 



?3 



-ci e 



i(Ci+^l 1) ) 



cos^^tanh^i^ + ) + isin^^ 



where 



(2) 

V«i y 



a) 



fclfl - |Ci| 2 COS 2 0i 

^(laf^ + laf'l 2 ) 



(ii) Soliton 2 (r? 2 i? - 0,ViR -> oo): 



V? 2 y 



/a 2 -\ 



V^y 



A^sech ( //,/; • ^ „ ' j ) 



93 -cie 



(fcj + i&i) 



cos^^tanh ^77 2j r + 



(i?3 — i?i 



+ isin02 1 ' ) 



where 



2 I Vl-«i«2 



V 



Q. 



(2)/ l-(gt 2) /4 2) ) Kl 



b) After collision (limit z — > oo): 



/ 



(i) Soliton 1 (77U2 ~ 0,r] 2R oo): 



(fc 2 - fcQjfci + V^ 2 r ~ IciPcos 2 ^^ 



g 3 ~q e 



i(Cl+ ^)) (fc 2 ~ ^i) 



cos0^tanh yrjiR + 



+ isin^ 



(i) 
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where 



/ (i) f i-c4"A.i")~A \ 



V 



v/1 — K1K2 



(ii) Soliton 2 (^ 2j r ~ 0,^ -> -00) 
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(h - k 2 )\h + k* 2 



k\ R — |ci| 2 cos 2 0^ 1 ' ) 



(fci + k* 2 )\kx - k 2 \ k 
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zsmi 
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where 



/ n 0) \ 
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(2) 
V «2 / 



IciPcos 2 ^ 1 ^ 



(2) 12 



(B4c) 



In the above equations <p^' and are defined in Eqs. ( lAldl) and (1A2dl) . respectively. 
Note that A l + ^ A l f, and also in general the intensities |^ + | 2 ^ \A l ~\ 2 , j,l = 1,2, except 

(1) (2) 

when = ^h) ■ The explicit relation between the intensities before and after collision is 

C*2 a 2 

given by Eq. (|22|) in the text. The phase-shifts $1 and $2 suffered by solitons Si and £2, 
respectively are given in Eq. (23a). 



2. 1-bright— 2-dark soliton solution 



Considering the 1-bright— 2-dark two soliton solution (17), the analysis in the asymptotic 
limits can be performed as follows (with km, k 2 R > and ku>k 2 i)- 
a) Before collision (limit z — > —00): 
(i) Soliton 1 (r] 1R ~ 0,^ -> -00): 



Qi 



Al-kme irilI sech[r]m + 



-c e ife+ *i J)) 



Ri 



( os0^tanh( rjm + ) + ismcp^ 



J = 1,2, 
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(B5b) 
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where 
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(ii) Soliton 2 {rj 2R ~ 0, r] 1R -» oo) 
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COS< 



i0) 



(i? 3 — i?i 



+ isin02 
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_i / fai - bj 



Qi = {fa - fa 



Q 2 = {fa + k* 2 



1 + 



?2R 



+ 



{fa - ibi){k 2 - ibi) {fa - ib 2 ){k 2 - ib 2 ) 



+ 



lr 2 l 

l c 2l 



(A;i-i6i)(A5 + «6i) {fa - ib 2 ){k* 2 
Here A 2+ is defined in Eq. ( 1B8cj) given below, 
b) After collision (limit z — > 00): 
(i) Soliton 1 [r) 1R ~ 0,r) 2R -> 00): 



Z&9 



1/2 



1/2 



^ + /c 1R e ir ' 1/ sech(^ + 



(i? 3 — i? 2 



- Cj e 



{k 2 + ibj) 



cos^ftanhl r] 1R H 



(7) 

+ isin0^ 
3 = 1,2," 



where 



41+ _ ,1- / QiQ^ 

(ii) Soliton 2 (^ ~ 0,77^ -> -00): 

Qi - Al + k 2R e^ sech (r) 2R + f) 



~ _ c . e i(o+4 J) ) 



cos02 t ann ( ^2_r + ) + ism<f>2 



U) 



J = 1,2, 



where 



a 



(!)• 1 



(1)* 

a 



j—\ k 2R - ( (ciPcos 2 ^^ + |c 2 | 2 cos 

^2_R V ^ 
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Note that in the above expressions \A l ^\ = \A[ \,l = 1,2, and thereby confirming that no 
intensity redistribution occurs in the case of bright soliton, and also for dark solitons [see 
Eqs. (B5b), (B6b), (B7b) and (B8b)]. 
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|cj=0.3 



|cj=0.1 



|cj=0.9 




FIG. 1: Intensity plots of one-soliton of the mixed CNLS equations with N = 2 for different 
values of the background parameter c\ for a fixed value of z. Note that the intensity of the bright 
soliton increases as the depth of the dark soliton decreases. The other parameters are chosen as 
ki = 1 + i,a ( l ) = l,6i = 0.2. 
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FIG. 2: Intensity plots of one-soliton of the mixed CNLS equations with N = 3 for different 
values of the a parameters for a fixed value of z. (a) = 0.2 + 0.01i,a^ = 0.5 + 0.05a, (b) 
= 1 + i, = 2 + i and (c) = 13 — 13i, ot^ = 17 + 0.31 The role of a parameters are seen 
both in the intensity and phase of the bright soliton while it affects the phase (central position) of 
the dark soliton. The parameters ci, &i, and k\ are chosen as |ci| = 0.56, b\ = 0.2, k\ = 1 + i. 



31 





FIG. 3: Elastic collision of (bright-dark) two-solitons in the mixed CNLS system for the N = 2 
case. The parameters are chosen as given in the text. 
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FIG. 4: Intensity profiles of two colliding bright solitons of the mixed CNLS equations with N = 2, 
before (z=-4) and after (z=4) collision: (a) in the absence of dark component (ci = in Eq. (11)); 
(b) in the presence of dark component (c\ ^ in Eq. (11)). The figure is plotted for special choice 
of parameters (as given in the text) with |ci| = 0.56. Note the elastic nature of the collision. 
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FIG. 5: Shape changing collision of two-solitons in the mixed CNLS system for the N = 3 case. 
The parameters are as given in the text. 
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FIG. 6: Intensity profiles showing the collision scenario of two bright solitons of the mixed CNLS 

a (i) a (2) 

equations with N = 3 in the absence of dark component, (ci = 0) with -4^ ^ -fey. (a) z=-4 and 
(b) z=4, given by the specific choice of parameters as given in the text. 
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FIG. 7: Intensity profiles showing the collision scenario of two bright solitons of the mixed CNLS 

a m a (2) 

equations with N = 3 in the presence of dark component {c\ / 0) with-fjy / -fey. (a) z=-4 and 
(b) z=4, given by the specific choice of parameters as given in the text. 
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FIG. 8: Intensity profiles showing the collision scenario of two bright solitons of the mixed CNLS 

a W Q (2) 

equations with N = 3 in the absence of dark component (ci = 0), for -4ry = -W. (a) z=-4 and (b) 
z=4. The parameters are chosen as given in the text. 
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FIG. 9: Intensity profiles showing the collision scenario of two bright solitons of the mixed CNLS 

a (l) (2) 

equations with N = 3 in the presence of dark component (ci / 0), for -4jj = -fey. (a) z=-4 and 
(b) z=4. The parameters are chosen as in the text. 
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FIG. 10: Intensity profiles showing the propagation of two dark solitons in mixed CNLS system, 
before (z=-5) and after (z=5) collision: (a) for the special choice of parameters (as given in the text) 
with -^y / — ^y; (b) for the special choice of parameters (as given in the text) with — = -^y. 

0^2 ^~^2 ^~"^2 
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